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Abstract 

Spectral distribution theory, which can be used to compare microscopic interactions 
over a broad range of nuclei, is applied in an analysis of two modern effective in- 
teractions based on the realistic CD-Bonn potential for OhQ no-core shell model 
calculations in the fp shell, as well as in a comparison of these with the realis- 
tic shell-model GXPF1 interaction. In particular, we explore the ability of these 
interaction to account for the development of isovector pairing correlations and col- 
lective rotational motion in the fp shell. Our findings expose the similarities of 
these two-body interactions, especially as this relates to their pairing and rotational 
characteristics. Further, the GXPF1 interaction is used to determine the strength 
parameter of a quadrupole term that can be used to augment an isovector-pairing 
model interaction with Sp(4) dynamical symmetry, which in turn is shown to yield 
reasonable agreement with the low-lying energy spectra of 58 Ni and 58 Cu. 
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1 Introduction 



Realistic NN potentials, whether derived from meson exchange theory (e.g., [I]) 
or chiral effective field theory (e.g., [2]), and their effective interaction deriva- 
tives, provide no a priori indication regarding how well they may or may not 
reproduce prominent features of nuclei, such as pairing gaps in nuclear energy 
spectra or enhanced electric quadrupole transitions in collective rotational 
bands, until actually employed in shell-model calculations. While such calcu- 
lations are often laborious and model dependent, a simple and straightforward 
evaluation of an interaction can be made using spectral distribution theory 
[314"] . Indeed, spectral distribution methods can yield a deeper understanding 
of the nature of an interaction and above all, its role in the development of col- 
lective and correlated many-nucleon motion [51loTITll8ll{?lllUfllfl2| . In particular, 
spectral distribution theory can be used to show through correlation coeffi- 
cient measures the similarity of interactions. Such analyses are independent 
of the averages of the interactions and yield an overall comparison across a 
broad domain of nuclei beyond what can be achieved by overlaps of nuclear 
states or detailed comparisons of two-body interaction matrix elements. 

In this paper we examine three modern /jo-shell interactions, specifically, two 
interactions denoted as "CD-Bonn" and CD-Bonn+3terms [13] based on the 
CD-Bonn realistic potential pQ as well as GXPF1 [H]. The GXPF1 effec- 
tive interaction is obtained from a realistic G-matrix interaction based on the 
Bonn-C potential [13] by adding empirical corrections determined through 
systematic fitting to experimental energy data in the fp shell. The CD-Bonn 
potential is a charge-dependent one-boson-exchange nucleon-nucleon (NN) 
interaction that is one of the most accurate in reproducing the available 
proton-proton and neutron-proton scattering data. Specifically, we use two- 
body matrix elements of an effective "CD-Bonn" interaction derived from the 
CD-Bonn potential for OhQ no-core shell model (NCSM) calculations in the 
fp shell. In addition, the CD-Bonn+3terms interaction introduces phenomeno- 
logical isospin-dependent central terms plus a tensor force with strengths and 
ranges determined in OhQ NCSM calculations to achieve an improved descrip- 
tion of the A = 48 Ca, Sc and Ti isobars. In the regard, we use spectral 
distribution theory to provide an assessment of differences between the novel 
CD-Bonn+3terms interaction and "CD-Bonn" as well as a comparison of these 
interactions with GXPF1, which has been shown to reproduce nuclear energy 
spectra throughout the fp shell [T4"] . 

The likely success of the three interactions for reproducing pairing and rota- 
tional spectral features is examined in a comparison to a pairing-plus-quadrupole 
model interaction, which combines a Sp(4) dynamically symmetric model in- 
teraction [ToTfTT] for description of like-particle and proton-neutron (isovec- 
tor) pairing correlations with a SU(3) symmetric term that governs a shape- 
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determined dynamics. If the model and effective interactions are strongly cor- 
related, then the latter will reflect the characteristic properties of the simpler 
model Hamiltonian, which in turn may be used as a good approximation. 

The present study, which is complementary to a similar If 7/2 analysis |18j . 
focuses on the upper fp-shell domain, which includes neutron-deficient and 
N m Z nuclei along the nucleosynthesis rp-path and unstable nuclei currently 
explored in radioactive beam experiments [T9"f2T)] . The analysis of the upper 
fp results reveals overall properties of the interactions under consideration 
different from those observed in the I/7/2 orbit. 

Several detailed reviews of the nuclear shell model and its applications have 
been published recently [2T|22|23j that delve into related key physics issues 
that are explored in the present work. However, the spectral distribution anal- 
ysis provided here is novel and sheds considerable light on the features of new 
fp-shell interactions, some of which have been developed since those reviews 
were completed. 



2 Theoretical Framework 

The theory of spectral distributions is an excellent approach for studying mi- 
croscopic interactions [4"|24II25] and continues to be a powerful concept with re- 
cent applications in quantum chaos, nuclear reactions and nuclear astrophysics 
including studies on nuclear level densities, transition strength densities, and 
parity/time-reversal violation (for example, see [26l[27p 8,29,30,31,32,33,3 4l[35] ). 
The significance of the method is related to the fact that low-order energy 
moments over a certain domain of single-particle states, such as the energy 
centroid of an interaction (its average expectation value) and the deviation 
from that average, yield valuable information about the interaction that is 
of fundamental importance [7|11|25II36|37 .38 39 4 0|^T] without the need for 
carrying out large- dimensional matrix diagonalization and with little to no 
limitations due to the dimensionality of the vector space. Note that if one 
were to include higher-order energy moments, one would obtain more detailed 
results that, in principle, should eventually reproduce those of conventional 
microscopic calculations. 

Spectral distribution theory (see Appendix for basic mathematical definitions 
and notation introduction) combines important features, the most significant 
of which are as follows: 

(1) The theory provides a precise measure, namely, the correlation coefficient, 
for the overall similarity of two interactions. Literally the correlation co- 
efficient is a measure of the extent to which two interactions "look like" 
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(are correlated with) one another. In this respect, correlation coefficients 
can be used to extract information how well pairing/rotational features 
are developed in interactions, which may differ substantially from an in- 
dividual comparison of pairing/quadrupole interaction strengths [T8] . 

(2) It gives an exact prescription for identifying the pure zero- (centroid), one- 
and two-body parts of an interaction under a given space partitioning. 
Therefore, major properties follow: 

(a) The correlation coefficients are independent of the interaction cen- 
troids. (A direct comparison of two-body matrix elements provided 
by AW potentials may be misleading, especially when the averages 
of the interactions differ considerably.) 

(b) The pure one-body part of an interaction, the so-called induced 
single-particle energies (Eq. [T4|) . is naturally identified in the frame- 
work of spectral distribution theory and is indeed the average monopole 
interaction (compare to [33]). As such it influences the evolution of 
the shell structure, shell gaps and binding energies |4"2"] . 

(c) The pure two-body part is essential for studies of detailed property- 
defining two-body interactions beyond strong mean-field effects. 

(3) The correlation coefficient concept can be propagated straightforwardly 
beyond the defining two-nucleon system to derivative systems with larger 
numbers of nucleons [I] and higher values of isospin [6]. This, in addi- 
tion to the two-nucleon information provided by alternative approaches 
(e.g.,[43j), yields valuable overall information, without a need for carrying 
out extensive shell-model calculations, about the universal properties of 
a two-body interaction in shaping many-particle nuclear systems. 

Group theory underpins spectral distribution theory [3l4ToTI£ll4"4] . The model 
space is partitioned according to particular group symmetries and each subse- 
quent subgroup partitioning yields finer and more detailed spectral estimates. 
Specifically, for n particles distributed over 4S1 single-particle states, the spec- 
tral distribution averaged over all n-particle states associated with the U(4f2) 
group structure is called a scalar distribution (denoted by "n" in the formu- 
lae) and the spectral distribution averaged over the ensemble of all n-particle 
states of isospin T associated with U(2f2) (g) U(2)t is called an isospin- scalar 
distribution (denoted by "n, T"). 

For a spectral distribution a (a is n or n, T), the correlation coefficient between 
two Hamiltonian operators, H and H', is defined as 




((H* - (H*) a )(H' - (H') a )) 



a 



a 



(1) 



(r H a H i 



(T H o- H > 



where the "width" of H is the positive square root of the variance, 



K) 2 = <(#-<#n 2 > 



(2) 
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and the steps for computing these quantities are outlined in the Appendix. 
The average values, (0) a , related to the trace of an operator O divided by the 
dimensionality of the space, are given in terms of the ensemble considered. In 
the (isospin-)scalar case, the correlation will be denoted by ( n (C"' T ) or simply 
C (C T ) f° r n = 2. The significance of a positive correlation coefficient is given 
by Cohen [15] and later revised to the following table: 

Table 1 

Interpretation of a correlation coefficient. 

trivial small medium large very large nearly perfect perfect 

0.00-0.09 0.10-0.29 0.30-0.49 0.50-0.69 0.70-0.89 0.90-0.99 1.00 



From a geometrical perspective, in spectral distribution theory every interac- 
tion is associated with a vector and the correlation coefficient ( (Eq. [I]) defines 
the angle (via a normalized scalar product) between two vectors of length a 
(Eq. [2]). Hence, Ch,h> gives the normalized projection of H onto the H' inter- 
action (or H' onto H). In addition, ((h,h>) 2 gives the percentage of H that 
reflects the characteristic properties of the H' interaction. 

The pairing and rotational characteristics of an interaction can be probed 
through its projection onto a model microscopic Hamiltonian that describes 
isovector pairing correlations and collective rotational excitations. While the 
latter possess a clear SU(3) symmetry [19] within the framework of the har- 
monic oscillator shell model, the former respect a Sp(4) dynamical symmetry 
[4"5|47|48| . Specifically, we employ the pairing-plus-quadrupole model interac- 
tion 

H M = # s p(4) + Hq(2), H q = - '-jQ ■ Q, (3) 

where H sp ^ is an isoscalar Sp(4)-dynamically symmetric interaction [17] for a 
system of n valence nucleons (an eigenvalue of N) in a 4f2-dimensional space, 

# SP (4) = -G £ A\K - f (f 2 - f ) - - eN, (4) 

i=-l 

with two-body antisymmetric JT-coupled matrix elements (Appendix Eq. [6]) 
for {r < (s, t); t < u} orbits, 

W r J L = -Go^5 JO S T1 5 rs 6 t u - {-E [(-) T + \\ + C}S rt S su , (5) 

where Go = G + j, E = (J^ + y), G, F, E, D and C are interaction strength 
parameters and e > is the Fermi level energy (see Table I in Ref.fTT] 
for parameter estimates). The sp(4) algebraic structure is exactly the one 
needed in nuclear isobaric analog + states to describe proton-proton (pp), 
neutron-neutron (nn) and proton-neutron (jpri) isovector pairing correlations 
(accounted by the pair annihilation (creation) operators A+\ _ 1 ) and isospin 
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symmetry. The latter is reflected by the isospin operator T 2 and related to 
a J-independent isoscalar (T = 0) pn force. The most general model inter- 
action with Sp(4) dynamical symmetry [I7J has been found to provide for 
a reasonable microscopic description of the pairing-governed isobaric analog 
+ states in light and medium mass nuclei and to account quite well for 
the observed detailed structure beyond mean-field effects such as the N = Z 
anomalies, isovector pairing gaps and staggering effects |16j . In addition, the 
i/ sp (4) interaction (j3J) has been shown to strongly correlate, especially when the 
quadrupole term is introduced, with the "CD-Bonn", CD-Bonn+3terms and 
GXPF1 interactions in the I/7/2 orbit [T8] . 

The Hq (2)-term in (Eq. [3]) is the part of the pure two-body Hq{2) quadrupole- 
quadrupole interaction that, in the vector algebra terminology, is orthogonal to 
the pure two-body Sp(4) Hamiltonian [7J. This is because the Sp(4) interaction 
contains a part of the quadrupole-quadrupole interaction that is not negligible 
as revealed by the correlation between Hq and H^y Namely, in the scalar 
case it is 15% (I/7/2), 29% (l/ 5/2 ) and 29% (2p 1/2 2p 3/2 ), and for the T=l part 
of the interactions, it is 34% (l/ 7 / 2 ), 58% (l/ 5/2 ) and 58% (2p 1/2 2p 3/2 ). 

Such a Hamiltonian (Eq. [3]) does not affect the centroid of H s pU) because 
Hq(2) is traceless. In this way this collective interaction preserves the shell 
structure that is built into H sp ^ and established by a harmonic oscillator 
potential and as a result is favored in many studies [TJH]. 



3 Results and Discussions 

The similarity of the "CD-Bonn", CD-Bonn+3terms and GXPF1 interac- 
tions, which will be denoted as H , and their pairing/rotational characteristics 
can be tracked in many-nucleon systems [TU] through the propagation formu- 
lae (Appendix Eqs. [121 [16]) - The latter determine how the averages extracted 
from the two-nucleon matrix elements get carried forward into many-nucleon 
systems. This propagation of information is model-independent. 

The present investigation focuses on the upper fp-she\\ domain and is com- 
plementary to a similar I/7/2 analysis [18] (a few results from that study are 
presented in Section [37T1 for completeness). Such a partitioning of the fp oscil- 
lator shell follows naturally from a splitting of these two regions by a strong 
spin-orbit interaction. We examine the H (2) pure two-body part of the fp- 
shell interactions and how it is correlated to the model interaction (Eq. [3]). 
The latter, in addition to its centroid, is pure two-body in the upper fp model 
space because of the assumption for H sp ^ of constant Fermi level energy and 
fixed interaction strengths throughout the entire region. The significance of 
the correlation coefficients between pure two-body interactions [25] reflect the 
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fact that nuclear states, their collective properties and configuration mixing, 
are solely shaped by the pure two-body part of an interaction, while the one- 
body part, albeit of a considerable significance, trivially reorders the states in 
the nuclear energy spectrum. In addition, such analyses are free of the one- 
body influence including induced single-particle energies, which are related to 
the monopole interaction [4~3f42|H4| . and external single-particle energies. The 
latter are introduced when a core is assumed, as for the Ohfl 40 Ca-core shell 
model using the GXPF1 interaction. For the OhQ NCSM calculations with 
"CD-Bonn" or CD-Bonn+3terms, the two-body matrix elements specifying 
the particle-core interactions supplant the role of external single-particle en- 
ergies. These additional two-body matrix elements together with the external 
single-particle energies for GXPF1 are not included in the present analyses. 



In our study, we vary only x, the quadrupole strength parameter in (Eq. [3]), 
to find its optimal value (which is an exact solution) by maximizing the ( 
correlation coefficient [5U] between the model Hm interaction and the pure 
two-body part Hq(2) of each of the effective interactions under consideration. 
We do not alter the parameters of the Sp(4) model, which have already been 
shown in an appropriate domain of states to be valid for reproducing various 
quantities (such as binding energies and pairing gaps) and are in agreement 
with estimates available in literature [TolfTT] . 

In both scalar (Table [2]) and isospin-scalar (Fig. [T]) distributions, the Ch^(2) h m 
correlation with the pairing+quadrupole Hm interaction is stronger for GXPF1 
compared to "CD-Bonn" and CD-Bonn+3terms. Hence, other types of inter- 
actions that do not correlate with the pairing and quadrupole-quadrupole 
interactions (of fixed strength throughout the upper fp shell) comprise a rela- 
tively small part of the pure two-body GXPF1 interaction. They are weakest 
for the T = n/2 group of states (Fig. [T]) for all the three interactions. In 
addition, the results reveal that the symplectic Sp(4) dynamical symmetry 
of Hm, especially when compared to both CD-Bonn interactions and in the 
highest-isospin (T = n/2) group of states, is only slightly broken by Hq(2) 
(Tabled fifth row and Fig. [T]). The Sp(4) symmetry breaking is related to the 
correlation coefficient of Hq(2) with Hm or equivalently to the ratio of their 
norms [51], where {Ch u ,hM2)Y = 1 - (Ch u ,h. pW ) 2 = 1 - #s P (4) (Table [2]). 



Two correlation coefficients, which in the present study are independent of 
any interaction strength parameters, are of particular interest. Specifically, 
the isospin-scalar space partitioning is where the ability of an interaction to 
form correlated pairs and hence reproduce prominent pairing gaps is detected 
via (,Ho{2) h (4) • The capability of an interaction to describe rotational collec- 
tive motion, and hence to reproduce rotational bands and enhanced electric 
quadrupole transitions, can be detected via its correlation to the full Hq(2) 
quadrupole-quadrupole two-body interaction, Ch^(1) h q {2)- 
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The results show a large correlation with isovector pairing of the T = 1 part of 
the /p-shell interactions under considerations (n = 2), especially for GXPF1, 
and a good tendency towards development of pairing correlations in the T = 
n/2 states (Fig. [I]). In both scalar and isosp in-scalar cases (Tabled fourth row 
and Fig. [2]), the rotational features are more fully developed for GXPF1 and 
less for CD-Bonn+3terms and "CD-Bonn" (Fig. [2]) different from the outcome 
in the I/7/2 region especially for the T = 1 part of "CD-Bonn". 

In short, the GXPF1 interaction is expected to reproduce spectral features like 
pairing gaps and rotational bands observed in the upper fp nuclei, while it is 
unlikely for both CD-Bonn interactions to fully reflect the rotational properties 
of these nuclei. In comparison, the CD-Bonn+3terms interaction in the I/7/2 
orbit exhibits well-developed pairing and rotational characteristics [18]. Such 
a difference in the behavior of CD-Bonn+3terms within both regions, 1/7/2 
and upper fp, may reflect the fact that this interaction was determined through 
a reproduction of the low-lying energy spectrum and binding energy of A = 48 
I/7/2 nuclei. 

The different extent to which the GXPF1 interaction compared to the "CD- 
Bonn" and CD-Bonn+3terms interactions reflects development of pairing cor- 
relations and collective rotational modes in the upper fp domain may be the 
reason why their pure two-body part do not correlate strongly as, for example, 
"CD-Bonn" and CD-Bonn+3terms do. Namely, in the (isospin-) scalar case 
the pure two-body correlations^ are 0.90 (0.88) between "CD-Bonn" and CD- 
Bonn+3terms and only 0.56 (0.37) between "CD-Bonn" and GXPF1 and 
0.53 (0.40) between CD-Bonn+3terms and GXPF1. In addition, the isospin- 
scalar correlation coefficients involving the significant induced pure one-body 
(monopole) contribution, Xj (Appendix Eq. [141) . differ between GXPF1 and 
the two "CD-Bonn" interactions. Their behavior, especially below mid-shell, 
reflects the similarity of the corresponding T = induced single-particle en- 
ergies and the opposite signs of the corresponding Ai^* and AfjT* . 

In summary, when compared to the original "CD-Bonn", Ohfl NCSM calcu- 
lations with CD-Bonn+3terms for upper fp nuclei are likely to achieve an 
improved description of many-body spectral phenomena, associated with pair 
formation (especially when T ^ n/2) and with rotational motion (especially 
in the highest-isospin states). While such results lie slightly closer to what is 
achieved by the GXPF1 interaction with two-body matrix elements directly 
adjusted to experimental /p-shell data, larger NSCM model spaces for both 
CD-Bonn interactions may be necessary for a better reproduction of pairing 
and rotational spectral features throughout the upper fp domain. 



2 In the isospin-scalar case, the correlations vary slightly with the particle number 
and isospin and average values are quoted. 
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3.1 Individual- Orbit Analysis 



One can further perform a partitioning of the /jo-space to single- j orbits, I/7/2, 
I/5/2, 2pi/ 2 and 2p 3 / 2 , to provide for more detailed spectral measures that may 
reflect important fine effects that are otherwise averaged out when the entire 
fp major shell is taken into account. We have already illustrated such an 
example and its significance by exploration of the 1/7/2 orbit [IS]. Individual 
orbit analyses render correlation coefficients that are free of the influence of 
any one-body interaction [by definition, (Appendix Eq. ITf9l[T4|) ]. However, due 
to the small model space, the 2pi/2 and 2£>3/2 orbits are combined and in their 
joint region only the pure two-body part of the interactions is considered. 

In the scalar distribution, a good portion, 53% to 98%, of the H M pair- 
ing+quadrupole model interaction is described solely by the inter- 
action (i?s P (4), Table [3]). The Hq /p-shell interactions exhibit a quite well- 
developed rotational character (Table [31 (h .h q ) except for "CD-Bonn" and 
CD-Bonn+3terms in the 2p!/22p3/ 2 region. Besides these cases, the H M model 
interaction, as revealed by (h ,h m i n Table El can be used as a very good ap- 
proximation to the /p-shell interactions within each of the domains considered. 



The (Jj o H rp ( 4) correlation coefficients (Table H]) show large J = isovector 
coherence within each single- j shell, particularly for the T = 1 part of CD- 
Bonn+3terms and GXPF1, which are expected to describe quite well phe- 
nomena of a pairing character, while for "CD-Bonn" other types of inter- 
action compete with pair formation. The latter are of the Hq{2) type for 
the I/7/2 and I/5/2 orbits, where the residual interactions are negligible. 
In short, the simple Sp(4) model interaction and especially its extended pair- 
ing+quadrupole Hm interaction can reproduce reasonably well the T = 1 
part of the three effective interactions under consideration within the orbits 
specified in Table EH 

Within individual orbits a very close similarity is observed between the effec- 
tive and model interactions as well as in nuclear systems with more than two 
nucleons. However, more prominent differences among the interactions appear 
in the multi-j upper fp domain especially concerning both "CD-Bonn" in- 
teractions. This may indicate that the inter-orbit interactions do not respect 
strongly the symmetries imposed in the model interaction. In addition, the 
interaction strengths may differ from one orbit to another. While they do not 
affect correlation coefficients in the singe- j cases, their relative strength is of 
a great importance for multi-j analysis. 
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3. 2 Energy Spectra for A = 58 Nuclei 



The results presented above show that the Hm model Hamiltonian (Eq. [3J can 
be used in the upper fp region as a quite good approximation of the pure two- 
body part of the GXPF1 interaction (Tableland Fig. [p. This implies that 
both interactions are expected to yield energy spectra of a similar pattern. 

As an illustrative example, we apply the simpler model interaction (Eq. [3]) to a 
nuclear system of two nucleons in the upper fp region without any parameter 
variation. Particularly, we assume a 56 Ni-core and that both nucleons in 58 Ni 
and 58 Cu occupy the upper fp orbits with reasonable probability [Hj . For a 
description of the low-lying structure of these nuclei, external single-particle 
energies are needed to rescale at the end the eigenvalues of the model Hamilto- 
nian. This is performed trivially due to the microscopic structure of the model 
eigenstates, which are constructed in terms of fermion creation operators. We 
adopt single-particle energies that are derived from the 57 Ni energy spectrum. 
To a very reasonable degree, these energies reflect the influence of the I/7/2 
orbit and the core mean-field contribution. 

The model Hamiltonian uses a x value of 0.027 that we obtained through a 
comparison of Hm to GXPF1 within a scalar distribution. The reason is that 
the pure two-body GXPF1 interaction holds the best correlation coefficient 
to the model interaction in the upper fp region (Table [2]). In addition, the 
energy spectra for 58 Ni and 58 Cu are found to be closely reproduced by shell 
model calculations with the GXPF1 interaction in the full fp shell |14j . 

We extend the model space to include the lgg/2 orbit as it intrudes in the 
upper fp domain. This is exactly the space where the Sp(4) model was applied 
and interaction strength parameters determined [17]. The results (Fig. [3]) show 
a very good reproduction of the low-lying T = 1 spectra in 58 Ni and 58 Cu, 
especially the lowest 2 + states for both nuclei and the first + (T = 1) state 
above the 58 Cu ground state. Both states are of particular significance. On 
the one hand, the energy difference between the lowest + and 2 + states is 
believed to be directly affected by the formation of correlated pairs in the 
lowest + state (ground state for even-even nuclei) and the pairing gap that 
occurs below the first excited 2 + state of a broken pair. On the other hand, 
the + (T = 1) to 1 + (T = 0) energy difference in N = Z odd-odd nuclei is 
associated with the close interplay of isovector (T=l, pairing correlations) and 
isoscalar (T=0) interactions between protons and neutrons in the same major 
shell. In addition, the T = spectrum of 58 Cu as predicted by the model 
interaction possesses the same pattern of the levels observed, namely, the 1 + 
ground state is followed by 3 + , 1 + , 2 + , 4 + , and 3 + . The 58 Cu T = spectrum 
appears narrower than the experimental data, which suggests that different 
quadrupole strengths for T = and T = 1 need to be used. 
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In short, we demonstrate that strong correlations typically yield very similar 
energy spectra and reproduction of the overall pattern of the energy levels 
without any adjustment of the interaction strength parameters. 



4 Conclusions 

In the present article, we have compared three modern /p-shell interactions, 
"CD-Bonn" , CD-Bonn+3terms, and GXPF1, based on realistic nucleon-nucleon 
potentials, in the upper fp region by means of the theory of spectral distribu- 
tions. We focused on the weaker but property-defining two-body part of the 
interactions and studied their pairing and rotational character in a comparison 
to a pairing+quadrupole model interaction, which includes Sp(4) dynamically 
symmetric isovector pairing correlations and a proton-neutron isoscalar force 
together with a quadrupole-quadrupole interaction for description of SU(3) 
dynamically symmetric collective rotational mode. 

The outcomes show that the pure two-body GXPF1 interaction, which is 
adjusted through systematic fitting to experimental energy data in the fp 
shell, demonstrates well developed collective rotational features and a ten- 
dency of its T = 1 part towards formation of correlated pairs. In addition, it 
strongly correlates with the pairing+quadrupole Hm model interaction, which 
is reflected in the quite good agreement between the experimental low-lying 
energy spectra of 58 Ni and 58 Cu and the theoretical prediction based on Hm 
in the l/5/22pi/22p3/2lg9/2 major shell. Individual-orbit analysis, including 
the I/7/2, I/5/2, 2pi/ 2 , and 2p 3 / 2 levels, shows considerably stronger correlation 
to Hm (up to 0.8 — 1.00) and a clear pairing or/and rotational character for 
all the three /p-shell interactions, which in addition correlate more strongly 
among themselves. In this respect, inter-orbit interactions may be a reason 
why both CD-Bonn interactions suitable for OhQ NCSM /p-shell calculations 
invoke small to medium pure two-body correlations in the upper fp domain. 
While these correlations for CD-Bonn+3terms appear slightly superior to the 
ones for "CD-Bonn" due to a phenomenological correction, it may be interest- 
ing to investigate how they vary for effective interactions based on the realistic 
CD-Bonn potential when higher- configurations are included in the no-core 
shell model analysis. 

In summary, based on these results, spectral distribution theory appears to be 
a good framework for uncovering fundamental properties of realistic interac- 
tions and their effective interaction derivatives in many-nucleon systems. We 
find varying degrees of respect for selected underlying symmetries. As some of 
these symmetries have been demonstrated to be important for certain spectral 
features, we have a tool for rapidly assessing the likely success of these inter- 
actions for reproducing those spectral features. For example, it is unlikely that 
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the CD-Bonn+3terms interaction will provide a fully satisfactory description 
of the rotational properties of nuclei in the upper fp shell. Given that this 
interaction was determined only with A = 48 nuclear spectra and binding en- 
ergies [13] , future efforts at expanding the region of its validity in the no-core 
shell model should benefit from the analysis provided here. 
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Appendix 

The theory of spectral distributions (or statistical spectroscopy) is well doc- 
umented in the literature [3l4l3"TfoTl8"] and is accompanied by computational 
codes [811512] for evaluating various measures. The purpose of this appendix is 
to specify the notation and ensure that our definitions of the summations and 
numerical factors that enter into such measures are clearly understood. 

In standard second quantized form, a one- and two-body interaction Hamilto- 
nian is given in terms of fermion creation a^nm^ = c J m (i/2W an( ^ annihilation 
%_ m (i/2)-CT = (—iy~ m+1 ^ 2 ~ a Cjm(i/2)cT tensors, which create or annihilate a par- 
ticle of type a = ±1/2 (proton/neutron) in a state of total angular momentum 
j (half integer) with projection m in a finite space 2Q = £j(2j + 1), 

^=-Efr{4®«J (00) 

r<s 

-\ E ^(1 + <W(1 + ttu)[T]Wrstu{Wr ® a\} T ® K ® a M } r }(°°\(6) 

rstuF 

where the labels are r = {j r , r r = |}, [r] = 2(2j r + 1), and [r] = (2J+ 1)(2T + 
1). In (Eq. ED, e rs is the (external) single-particle energy (hereafter we consider 
no angular momentum degeneracy for two different radial quantum numbers, 
e rs = e r 5 rs ) and W^ u =< rsJTMT \H\tuJTMT > is the two-body antisym- 
metric matrix element in the JT-coupled scheme [W^ tu = — (— ) r+s ~ r W^ rtu = 
_^_y+u-vy^v^ _ ^_y-+s-t-uyyT^ W^J. For an isospin nonconserving 

two-body interaction of isospin rank T, the coupling of fermion operators is 
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e r = e r - e = e r - — J2 £ sNs, (7) 



as follows, {{aj ® cl\} jt ® {a< <8> a«} JT } < - 0T ' ) , with H / vS l JT matrix elements. 

For a major shell that consists of several s orbits, each of degeneracy J\f s 
(M = J2s -A/"s): the (traceless) external single-particle energy of the r th orbit is 
obtained as 

1 

Af , 

where the average (external) single-particle energy is e = jj J2 S £sA/" s . 

Scalar Distribution. For a two-particle system, the monopole moment (cen- 
troid), which is the average expectation value of the two-body interaction, is 
defined in the scalar 

w - 1 V" rriM/ r - £™,r[r]W" r r 8rB (l + 6rs) , . 



where M — 4fi = 2E r (2jr + 1), the T-sum goes over all possible (J, T) for 
given r, s, and (^Q = X) r < Sj r[r]. The traceless induced single-particle energy 
is constructed by contraction of the two-body interaction into an effective 
one-body operator under the particular group structure, 

^ = T?E [JT]W r f rs (l + 8 rs ) 4E IJT]W£ U (1 + 5 tu ). (9) 

JVr s,JT yV tu,JT 

For a system with one hole in the r th orbit, A r corresponds to the energy of 
a single particle as contributed by the interaction with the valence particles 
above the core. In turn, the traceless pure two-body interaction is defined as 

w£M = Kil - (w c + ^^)5 rt 5 su . (io) 



In order to calculate energy moments and their propagation for higher n (and 
T) values, each interaction H (consisting of one(fc = 1)- and two (A; = 2)-body 
parts) needs to be expressed as a linear combination of terms of definite par- 
ticle rank (irreducible tensors 'Hkiy) of rank u = 0,1,2), that IS dbS db collection 
of pure zero-, one- and two-body interactions. For n particles, the Hamiltonian 
can be rendered (the sum J2* goes over r < s, t < u and Y = ( J, T)), 



H = nH 1 (0) + ™ Wo) + Hx{l) + (n - 1)W 2 (1) + H 2 {2) 



-ne 



2) W ° ~ IJ^dr + ^fA^'M ® O-}""" (ID 
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for then the quantity that defines the correlation coefficient (Eq. [T]) is easily 
computed for different particle numbers n, 



• ~ m ~ n){Ar ~ n ~ l) VmU2)W*^2). (12) 



Af(N - l)(jV - 2) (A/" - 3) 



Isospin-Scalar Distribution. Analogously, the centroid is defined as 

2 



' 7WTFm £j ms - (13) 



where M = 2Q. The AjT traceless induced single-particle energy for orbit r and 
the W^ u {2) traceless pure two-body interaction [8] are defined as, 



% = w T,lJ]w r J s l(i + <U - EiJW/l(i + 5 tu ), (14) 

JX r s,J yV tu,J 

W r J L(2) = W r J s l - (Wj + j^l^T ^rtSsu- (15) 

In order to calculate the correlation coefficient ( n,T and the variance a n ' T , the 
following quantities are needed (the sum J2* goes over r < s, t < u and J), 



(H^H') n ' T ~ (HY' T (H') n ' T = Pi{T)Y,Ke r e' r + £?i(n, T, r)M r [e r \' T r + KKl 

r r,T 

r,{ri,T2} 

+ £p 2 (n, T, r) ^ ElJ]^ u (2)^ r J ; fa (2), (16) 

where r is or 1, and the set {ti, r 2 } is {0, 0}, {0, 1} or {1, 1}. The propagator 
functions are derived in [531161 to be, 



m _ n(Af+2)(Af-f)-2AfT(r+l) , , 

J - AT(AT-1)(AT+1) K 1 ') 



n (n T t\ 4ArT(T+l)(l-»)(l-(-l)-)+(Ar + 2)(Ar-|)[( 2 r+l)n(n+2(-ir)-4T(T+l)(-l)^] 

WV'i J ) T ) — 4JV(Af-l)(Ar+l)(Af+2(-l)-r) 
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p 1 (n,T,r 1 ,r 2 ) = 8( ^ 1)(Ar+1)(Ar 1 _ 2)(Ar+2( _ iri) {4ArT(T+ l)(n - l)x 
(A/"-2n+4)(l-(-l) ri ) + [(2n + l)n(n + 2(-l) ri ) - 4T(T + l)(-l) ri ]x 

71 

[(2r 2 + l)(n + 2(-l)-)(AA - -)± + T(T + l)(-l)-][Ar - 2(-l)-]} 



„ („ T nr - C\\ - h("+2)-4T(T+l)][(Ar-f)(Ar-|+l)~T(r+l)] 

V2\n, i,t - u; - aA/w-i) 



ft (n, T, r = 1) = WTTJ ^ I ^ {lr 2 (T + 1) 2 (3AT 2 - 7JV + 6) 

i Tl Tl 

+-n(n-2)(Af--)(Af--+l)(Af+l)(Af+2) 

Tl 

+\T{T + l)[(5jV- 3)(M + 2)n(- - Af)+J\f(Af - l)(A/"+ l)(A/" + 6)]}. 



For the Sp(4) interaction, the average two-body interaction is expressed in 

terms of the model parameters in the scalar distribution as, W c = — + 

1 2 J 

jj^ppj ~ C, and in the isospin-scalar case for a given isospin value as, Wj = 
— 5ti+E q [(— l) T +i]— C. The pure one-body part of the Sp(4) Hamiltonian 



2(J\f—\^ — w 5 emu. in ijii<3 iDwopiii-o^ciicii ^aoc iwi a given louopn vaiuc cio, yy c 

-^5 T1 +E [(-l) T +l 

is zero for a singe- j orbit (by definition) as well as for the l/s/22pi/2 2^3/2 1^9/2 
major shell (by assumption of constant interaction strengths). The pure two- 
body matrix elements, W^ u (2), and hence the correlation coefficients involv- 
ing ifgp(4), are then independent of the C parameter in the scalar case and of 
the C and E parameters in the isospin-scalar case. 
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— •— (CD-Bonn.Sp(4)) 
— (CD-Bonn. H g ± (2)) 
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Fig. 1. Correlation coefficients of the pure two-body "CD-Bonn" (blue), 
CD-Bonn+3terms (red) and GXPF1 (green) interactions with -ff sp (4) (filled sym- 
bols), Hq(2) (transparent symbols) and Hm (empty symbols) in the upper fp shell 
for the isospin-scalar distribution. For each valence particle number, n, the isospin 
T varies as §,§ — 1, . . . , 0(|). The figures are symmetric with respect to the sign of 
n-2VL (n = 6). 



Table 2 

Scalar distribution correlation coefficients for many-nucleon systems of the Hq(2) 
pure two-body part of the "CD-Bonn", CD-Bonn+3terms and GXPF1 interactions 
with Hm (Eq.[3|), H sp ^ (Eq.[3J), Hq (2), and with the pure two-body full quadrupole- 
quadrupole interaction, Hq{2). The i?sp(4) = (CH M ,H ep ( 4) ) 2 quantity gives the part, 
in %, of Hm that is Sp(4) symmetric. 



'CD-Bonn" CD-Bonn+3terms GXPF1 



Ch {2),h m 

(h {2),H q (2) 



R 



Sp(4) 



0.57 
0.55 
0.14 
0.28 
93.1% 



0.54 
0.50 
0.20 
0.33 
85.7% 



0.83 
0.65 
0.51 
0.67 
61.3% 
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Fig. 2. Comparison between the orthogonal Hq{2) and the full two-body Hq(2) 
quadrupole-quadrupole interactions in their correlation to the pure two-body part 
of the /p-shell interactions, "CD-Bonn" (blue squares), CD-Bonn+3terms (red 
diamonds) and GXPF1 (green circles), in the isospin-scalar distribution. For each 
valence particle number, n, the isospin T varies as S, S — 1,...,0(±). The fi gure is 
symmetric with respect to the sign of n — 20 (0 = 6). 

Table 3 

Correlation coefficients for many-nucleon systems of the "CD-Bonn", CD- 
Bonn+3terms and GXPF1 interactions, Hq, with Hm (Eq. [3|) , -ff sp (4)(Eq. and 
the full quadrupole-quadrupole interaction, Hq. The -Rsp(4) = {Ch m ,h svW ) 2 quan- 
tity gives the part, in %, of Hm that is Sp(4) symmetric. 

Scalar Distribution 







!/7/2 


I/5/2 


2pi/ 2 2p 3 /2 a 


"CD-Bonn" 


(Hq,H m 


0.81 


0.76 


0.61 






0.66 


0.70 


0.58 




Ch ,H q 


0.69 


0.64 


0.17 




- R Sp(4) 


65.9% 


83.7% 


91.7% 


CD-Bonn+3terms 


(h ,H u 


0.87 


0.86 


0.52 




C-ffo,^sp(4) 


0.64 


0.67 


0.51 




Ch ,H q 


0.80 


0.82 


0.22 




- R Sp(4) 


53.4% 


60.6% 


98.3% 


GXPF1 


CHo,H M 


0.93 


0.84 


0.90 




C-ffo,^sp(4) 


0.76 


0.77 


0.70 




(h ,Hq 


0.78 


0.69 


0.85 




R Sp(4) 


67.9% 


84.7% 


61.7% 



pure two-body part of the interactions 
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Table 4 

Correlation coefficients for a two-nucleon system of the "CD-Bonn", CD- 
Bonn+3terms and GXPF1 interactions, Hq, with iJj^ (Eq. ED and H s pu\ (Eq. H]) . 
-Rs P (4) = (Ch m ,h sp(4) ) 2 gives the part, in %, of iT M that is Sp(4) symmetric. 



Isospin-scalar Distribution, T = 1 

I/7/2 I/5/2 2pi/ 2 2p3/2 



"CD-Bonn" 


xr=i 

^H ,H M 


0.95 


1.00 


0.59 




C T=1 

S #0:#sp(4) 


0.61 


0.57 


0.54 




- R Sp(4) 


41.5% 


32.3% 


84.2% 


CD-Bonn 


XT=1 


0.98 


1.00 


0.61 


+3terms 


C T=1 

S #0,# s p(4) 


0.85 


0.93 


0.59 




- R Sp(4) 


73.9% 


86.9% 


92.1% 


GXPF1 


r T=l 


0.96 


1.00 


0.94 




C T=1 

S #0,# s p(4) 


0.71 


0.86 


0.69 




-Rsp(4) 


54.5% 


74.0% 


53.5% 



pure two-body part of the interactions 
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Fig. 3. Theoretical ('th') low-lying energy spectra for 58 Ni (left, blue) and 58 Cu 
(right, red) compared to experiment ('exp', black). The theoretical calculations are 
performed in the 1/5/22^1/2 2^3/2 159/2 major shell with the Hm model interaction 
(Eq. [3] with x = 0.027) and with single-particle energies derived from 57 Ni experi- 
mental energy levels. 
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